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Abstract 



We consider the early cooling evolution of strongly magnetized strange stars in a CFL phase 
with high gap A > 100 MeV. We demonstrate how this model may explain main features of 
the gamma-ray burst phenomena and also yield a strong star kick. The mechanism is based on 
beaming of neutrino emission along the magnetic vortex lines. We show that for sufficiently high 
initial temperatures To ~ 30 to 60 MeV and surface magnetic fields B s ~ 10 15 to 10 17 G, the 
energy release within the narrow beam is up to 10 52 erg with a magnetic field dependent time scale 
Qh' between 10~ 2 s (for a smaller magnetic field) to 10 s. The above mechanism together with the 

parity violation of the neutrino-producing weak interaction processes in a magnetic field allow for 
the strange star kick. The higher the magnetic field the larger is the stars kick velocity. These 
velocities may cover the same range as observed pulsar kick velocities. 



X ■ 1 Introduction 



Compact stellar objects are subdivided into three groups: usual neutron stars (having only hadronic 
matter interiors and an extended crust), hybrid stars (with quark core, hadronic shell and a crust 
similar to that in neutron stars) and strange stars (hypothetical self-bound objects composed of u, d 
and s quark matter with only a thin crust if any). The possibility of strange stars has been extensively 
discussed in the literature, see [T] and Refs. therein. 

Gamma ray bursts (GRB) are among the most intriguing phenomena in the Universe, see |2] and 
Refs. therein. If the energy is emitted isotropically, an energy release of the order of 10 53 to 10 54 erg 
is needed to power the GRB. However, there is now compelling evidence that the gamma ray emission 
is not isotropic, but displays a jet-like geometry. In the case of beaming a smaller energy, of the order 
of 10 51 erg, would be sufficient to power the GRB |3 b . Although typically one may speak about a short 
GRB with a time scale t ~ 10~ 3 to 10" 1 s and a long GRB with a time scale t ~ 1 to 100 s, energies, 
durations, time scales involved in a burst, pulse shape structures, sub-burst numbers, etc, vary so much 
that it is hard to specify a typical GRB. 

Many scenarios for the GRB have been proposed: the "collapsar", or "hypernova" model linking 
the GRB with ultra-bright type Ibc supernovae (hypernovae) and the subsequent black hole formation 
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[3]; neutron star mergers [3], or accretion of matter onto a black hole; strange star collisions the 
model [7j assuming a large surface magnetic field up to 10 16 G of a millisecond pulsar produced in the 
collapse of a ~ 10 9 G white dwarf, with a powerful pulsar wind, as the source of the GRB; the model jS] 
suggesting a ee + plasma wind between heated neutron stars in close binary systems as consequence of 
the vv annihilation; the model JU] of a steadily accreting ~ 10 6 G white dwarf collapse to a millisecond 
pulsar with a ~ 10 17 G interior toroidal field, causing the GRB; an (isotropic) first order phase transition 
[TT] of a pure hadronic compact star to a strange star (see also |12j); asymmetric core combustion [12] 
in neutron stars due to the influence of the magnetic field generating an acceleration of the flame in the 
polar direction, etc. 

Pulsars are rotating neutron stars with rather high magnetic fields causing observable radio dipole 
signals. Most of the known pulsars are born in the neighborhood of the galactic plane and move away 
from it with natal kick velocities which are typically higher than those of their progenitors [T3|. This 
implies that the birth process of pulsars also produces their high velocities and thus cannot be entirely 
isotropic |15j . Up to now, the mechanisms driving this asymmetry are far from clear. There are several 
hypotheses, ranging from asymmetric supernova explosions , neutron star instabilities [T%1 ITU] and 
magneto- rotational effects [201 E] to the model of an electromagnetic or neutrino rocket It is also 
not clarified whether the distribution of pulsar kick velocities is bimodal with a low-velocity component 
of v < 100 km/s (20 % of the known objects) and a high- velocity component of v > 500 km/s (80%) 
as suggested by |22] or whether it can be explained by a one-component distribution JTH]. Most of 
the models are capable of explaining kick velocities of v ~ 100 km/s, but it is a nontrivial problem to 
explain the highest measured pulsar velocities around 1600 km/s. 

In this work we assume that a strongly magnetized strange star has been formed in the color 
superconducting color-flavor- locked (CFL) state with a large gap A(T), A(T = 0) > 100 MeV, cf. 
23J, 1 as the result of a phase transition. The latter could be caused by the accretion of the matter 
from a companion star, the neutron star angular momentum decrease owing to the gravitational and 
electromagnetic radiation after the supernova event had occurred, or during the collapse of a magnetized 
white dwarf to the neutron star state, the proto-neutron star collapse to the new stable state, or by some 
other reason. We also can deal with a hybrid star instead of the strange star. However, cf. there 
is experimental evidence that the GRB carries only a tiny baryon load of mass < 1O _4 M . Therefore 
the hadronic shell and the crust should be rather thin or there should be a special reason for a low 
baryon loading. We suggest to explain the beaming by the presence of a strong magnetic field. Thus the 
beaming and the low baryon loading stimulate us to conjecture about a magnetized strange star, which 
has a tiny hadron shell, if any, and only a thin crust. Varying the value of the surface magnetic field 
B s , we search for an optimal configuration to explain the GRB characteristics and to estimate a range 
of velocities of the star kicks. We use units h = c = 1. 

2 Vortices in magnetic field, neutrino beaming and GRB 

Vortex structures in the CFL phase of the strange star in the magnetic field. Let us further 
assume that the bulk of the star is in the CFL phase with no hadronic shell and a tiny crust generating 
a surface magnetic field which may achieve the values B s ~ 10 12 — 10 16 G. Then, one may expect 
by flux conservation that the inner magnetic field (at the crust-core interface) Bi U)S = B s (ni n /n s ) 2 ^ 3 
can reach the values <; 10 14 — 10 18 G, where n s and ni n are the corresponding densities at the surface 
and in the core, respectively. Below we will vary B in ^ s within the interval 10 15 < B in ^ s < 10 17 G. From 
the BCS relation the critical temperature of the superconductivity is T c = 0.57 A(T = 0) and the 
temperature dependent gap can be parameterized as [2S| A(T) ~ A(T = 0)[1 — (T/T c ) 3,4 ] a53 . For 
initial temperatures T ~ 30 — 60 MeV and the gap A(T = 0) <; 100 MeV we have T < T c . Thus the 

1 Our results are also relevant for the mCFL phase |23]. We only need all quarks to be gapped with large gaps 
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strange star is formed in the superconducting phase. Refs. J2IH 123 123 EH] have shown, that typical 
values of the Ginzburg-Landau parameter for color superconductors are k G l — A/£ ~ 1. The coherence 
length £cfl and penetration depth Acfl for CFL matter can be estimated in the weak coupling limit as 

*» * 0,(1^) (X - fin, W , ,( **) (^) (l - (1)")- fio , (1 ) 

where we used the above A(T) dependence. In our case the strong coupling constant a s = g 2 /47r ~ 1, 
and the quark chemical potential is fi q > 350 MeV. The critical Ginzburg-Landau parameter is k c gl ~ 
l/v2 distinguishes between type I (^gl < k gl) anc ^ type H ( k gl > k gl) superconductors. Then for 
T c > T 6 I_n ~ 15 -i- 20 MeV we deal with a type II superconductor. The latter entails the existence of 
a mixed phase (normal vortices embedded in superconducting matter) for T < T c in a broad interval 
of magnetic fields B cl < B < B c2 . The value B c2 is ~ B c k gl , with B c > 10 18 G for A > 100 MeV. 
Actually, even if B < B c \ = B c j kgl, the magnetic field, in spite of the Meissner effect, cannot be 
expelled from the star interior within a pulsar lifetime j30j , thus being concentrated in vortices aligned 
along the magnetic axis (feasibly being parallel to the rotation axis). Therefore, the mixed state exists 
for all B < B c2 . Making below only rough estimates we will assume the simplest Abrikosov vortex 
structure of the mixed phase. 

Due to magnetic flux conservation the number of vortices in the strange star interior is N vo = 
iiB injS R 2 /& q , where $ q = 6<3> with the magnetic flux quantum $ = 2 • 10~ 7 G cm 2 and the star radius 
R ~ 10 km I2E1 123 l2Ej • Superconductivity is expelled from the vortex interior, r < £, where the matter 
is in the state of a strongly magnetized quark-gluon plasma. The magnetic field decreases on the scale 
of the penetration depth Acfl as B ~ B exp(— t/Acfl), for r > Acfl, cf. [3TJ 1213 123 I2E] • Typical 
values of the magnetic field in the vortex center are Bq ~ ^q/i^^cFh) ~ 10 18 — 10 19 G, depending on 
Acfl- 

In hadronic matter and in metals, the pairing gap is rather small and the coherence length £ oc 
1/A(T) is larger than the Debye screening length. Therefore, charge neutrality is easily fulfilled for 
r < £. In the quark matter case we have opposite situation [321 EH]- For both, quarks and electrons, 
\ D = u\ C fl with typical values of v ~ 3 — 5. The presence of the Coulomb field diminishes the 
difference between the chemical potentials of quarks of different species since one does not need to fulfill 
the charge neutrality condition locally at distances r < but rather globally at r > A^. This implies 
that the cylindrical volume within distances r < from the vortex center is possibly filled by a 2SC+S 
phase plus a normal phase |34j . also characterized by a still rather high magnetic field B ~ 10 17 G. In 
case of the two-flavor superconducting (2SC) phase quarks of two colors, e.g., green and red, are paired 
with a large gap, whereas blue quarks are unpaired or paired with only a small gap < 1 MeV [23J. 
The 2SC+s phase emits neutrinos with a still high rate of about 1/3 of the emissivity of normal quark 
matter. In this phase at densities under consideration typical values of the electron fraction Y e = n e /nt,, 
where n e and n& are electron and baryon densities, are ~ 10~ 2 (for r ~ Ad)- The electron fraction 
in the CFL phase at T = vanishes [21], but at finite temperatures Y^ FL {T) ^ due to thermal ee + 
excitations. 

The typical distance between vortices is d = ($ g /5j„ jS ) 1//2 . The region near the magnetic vortex, 
where the magnetic field decreases by up to two orders of magnitude compared to the central value B , 
has the volume V vo ~ 2ti R{v\cfl) 2 ~ 10 21 fm 3 . The fraction of the total vortex volume N vo V vo /V (N vo 
is the total number of vortices, V is the star volume) varies from several -10~ 7 to several -10~ 3 for Bi n)S 
from ~ 10 13 G to 10 17 G, respectively. Therefore the condition d 3> Ad is safely fulfilled. 

Goldstone transport. The CFL (and mCFL) phase is also characterized by the presence of 
(almost) massless weakly interacting Goldstone excitations (like phonons). For T ^> ~ 2 — 3 

MeV of our interest the mean free path of Goldstones is much shorter than the star radius ~ 10 km, 
cf. jSn], and the typical time of their transport to the surface is very large. In the presence of a strong 
magnetic field generating vortices Goldstones cannot pass from CFL matter to the 2SC+S matter and 
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to the normal quark matter phase of the vortex core, where they cannot exist. Thus, the mean free 
path of Goldstone excitations Xq B in presence of vortices is still much shorter than in a homogeneous 
CFL strange star matter. For magnetic fields of our interest one finds Xq B ~ d 2 /£ « 6 x (10 2 — 10 6 ) fm. 
Thus Goldstones may efficiently equilibrate the temperature between the vortex interior and exterior 
regions. 

Neutrino radiation and beaming. The emissivity of the quark direct Urea (QDU) process in 
normal quark matter in absence of the magnetic field is very high 

e QDU ^ 2 x 1Q 26 Qs u y-1/3 rp6 QTg cm ~3 g -l ^ 

Correspondingly, the neutrino mean free path in normal quark matter becomes very short 

\® DU = 3 x 10 6 a; 1 u- 1 I7 1 / 3 T 9 " 2 cm. (3) 

For T <C ppe one Y e = pFe/PFb, where pp% = (^n 2 ^) 1 ^ 3 , i = e,p and u = nb/n is the compression. 
In superfluid matter the suppression factor of the emissivity is roughly [SZ| 

( DU ~ exp(-A(T)/T), T<T C , (4) 

and the enhancement factor of the mean free path is then Cdu- m ^he CFL phase pp e = but there 
exist thermal ee + excitations. Then pp e should be replaced by T and we may use above expressions 
(including approrpiate suppression and enhancement factors) with K e ~ T 3 / (37r 2 )nb. The emissivity 
remains virtually unaltered in the presence of a magnetic field when many Landau levels are occupied 
[38] and it is increased by 10 — 100 times if only the first Landau level is filled. However, the latter 
occurs only for B > 10 18 G. With the field values which we have assumed this condition is satisfied 
only in the region of the vortex cores, r < £, of a tiny volume. In the remaining part of the vortex 
volume V vo (£ < r < v\d, £ "C v\d) one can roughly neglect the influence of the magnetic field on 
eq. (J2J). Also multiplying (J2J) by a factor 1/3 < s < 1 we may estimate the emissivity of the 2SC+S 
+normal phase since due to the pairing in the 2SC phase the contribution of the red and green quarks 
is exponentially suppressed according (jlj). Also in the region between the vortices occupied by CFL 
phase with the volume V — N vo V vo , the emissivity of the direct Urea process is suppressed by (jH). The 
neutrino mean free path increased by this factor exceeds the radius of the star for temperatures below 
Topac ~ 20 MeV (for A(T = 0) ~ 150 MeV), whereas for the normal quark matter one gets T^° r c m < 
few MeV, cf. For higher gaps the value T^ C L still increases. 

For T > T££^ the heat transport is governed by the transport equation, CyT = kAT, Cy is the 
total specific heat of the matter and k = k v + kq is the total heat conductivity. Since Goldstones 
(phonons) are efficiently captured by vortices one has k„ > kg, and thus k ~ k v . A typical time scale 
for the energy transport to the surface is t tr ~ R 2 Cy/n. Since n u ~ Cy\ u , the energy transport time is 
t tr ~ R 2 /(\J). Using the above value \® DU for T ~ 20 - 40 MeV, we get t tr ~ 10~ 4 - lO^ 1 s (estimate 
is done for A(T = 0) ~ 150 MeV). For that time the star essentially losses its initial thermal energy. 
If initially (for T = 40 MeV, A > 100 MeV) the thermal energy is mainly concentrated in massless 
excitations, the energy falls by a factor 2 4 for T ~ T^ C L ~ 20 MeV. After the transport time t tr the 
CFL regions become transparent for neutrinos. 

For T < T p^ c L , the prominent difference between the neutrino mean free path in the 2SC+s+normal 
matter on the one hand and the CFL quark matter on the other hand leads to anisotropic neutrino 
emission via the direct (no rescatterings) QDU reaction from the star within a cone of the temperature- 
dependent opening angle Q v (T) ~ \ U (T)/R around the magnetic axis, where \ U (T) = \ U (T) V/(N vo V vo ). 
Thus after t > t tr the star begins to radiate the remaining energy ~ 10 52 erg within a narrow beam 
cone. 
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Cooling evolution. Numerical results. The cooling evolution for t > t tr within the opening 
angle can be described by inverting the solution of 



T(t) 

t = - J dT' C V (T')/L(T'), (5) 

rpCFL 

J opac 

where 

C V (T) = {[1 - N vo V vo /V] Cdu + sN vo V vo /V}C v (T) + C?(T), (6) 

Cy(T) ~ 10 39 M 2 / 3 (i?/10km) 3 Tg erg/K is the specific heat of the normal quark matter and Cy(T) ~ 3 ■ 
10 32 g ex (-R/10km) 3 T! eT g^ j s ^he contribution of (almost) massless Goldstone excitations and electrons 
and positrons (for the effective value of the degeneracy factor g ex ~ 11, we have assumed Vq — l/v 7 ^ 
for the velocity of Goldstones, the maximum Fermi velocity for relativistic fermi liquids. If vq were 
~ 1, we would get g ex ~ 3.). The luminosity is 

L{T) = [1 - cos e v {T)]L Q {T) ~ L (T)e 2 u (T)/2 (7) 

that takes into account the effect of the neutrino beaming. Here the isotropic luminosity is 

L (T) = N vo V vo se? DU (T, Y™) + (V - N vo V vo )e? DU (T, Y e CFL (T)) Cdu- (8) 

We take the electron fraction Y£° ~ 10 -2 for the 2SC+S + normal phase and a thermally produced 
electron fraction V e L (T) for the CFL phase. 

Here the Goldstones play an essential role since they rapidly (in a microscopic time) equilibrate the 
temperature between the vortices (of volume N vo V vo ) and their surrounding. In Eq. (JHJ) we neglected 
this very short time assuming an instantaneous heat transport from the hotter CFL regions of vortex 
exteriors to the regions of the vortex volume which then cool down by the direct neutrino radiation 
within the beaming angle. 

Results for the time evolution of neutrino luminosity, temperature and the neutrino beaming angle 
are shown in Fig. Q They demonstrate the effect of the cooling delay due to neutrino beaming in a 
strong magnetic field. For gaps ~ 100 — 200 MeV and B iritS ~ 10 15 G we get 9 U ~ 10 grad and the 
energy ~ 10 52 erg that remained in the star after the neutrino transport era has passed is radiated in 
the beam during several ~ 10~ 2 — 10 _1 s. Decreasing further the magnetic field one would get increase 
of the opening angle 9 V angle finally yielding inappropriate smearing of the beaming. For Sj H)S ~ 10 17 G 
we obtain a narrow beam, 9 V ~ 0.1 — 1 grad, and the initial star energy is radiated in the beam during 
several seconds. 

Neutrino energy conversion to photons, GRB. The neutrino/antineutrino collisions produce 
ee + pairs outside the star, which efficiently convert to photons. We use the result of jB] for the vv — > ee + 
conversion rate: 

E ee+ ~ 5 • 10 32 Tl 9 (R/ 10km) 3 erg s" 1 (9) 

to estimate that for a surface temperature T s ~ 10 — 30 MeV most of the initial thermal energy 
can be converted in ee + for an appropriate short time step t ~ 10 2 — 10~ 2 s. The advantage of the 
strange star model is that the strange star has a thin hadron shell and a tiny crust, if any. In this 
case we may assume that the surface temperature is of the same order of magnitude as the internal 
temperature. Note that in the presence of beaming already Ej ~ several • 10 51 erg produced on a time 
scale t ~ 10 -3 — 10 3 s could be sufficient to explain GRB. Thus, above estimates are very optimistic for 
the GRB model dealing with a magnetized quark star, if T s > 10 MeV. 

Moreover, ee + pairs are accelerated in the strong magnetic field of the strange star exterior, increasing 
the hard component of the X ray spectrum, in coincidence with experimental findings. 
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Figure 1: Time evolution of 
neutrino luminosity, temperature 
and neutrino beaming angle for 
diquark gaps A = 100 MeV 
(left panel) and A = 200 MeV 
(right panel), for different mag- 
netic fields B = Bi„ „. 



3 Strange star kicks from parity violation 

Range of kick velocities. Parity non-conservation in the weak interaction neutrino induced direct 
processes in presence of a magnetic field leads to a violation of reflection symmetry since the neutrino 
flux is a polar vector while the magnetic field is an axial one. The magnitude of this asymmetry has 
been estimated considering different reactions in neutron stars as the modified Urea reaction [39j, see 
also |40j : f3 decays on the pion condensate formation and breaking of nucleon pairs 02]; coher- 
ent neutrino electron scattering on nuclei J3H]; neutrino-polarized neutron elastic scattering; polarized 
electron capture [H]; the v — > vee + process jlH]; the reaction vv — > ee + ; etc. Some of these reactions 
produce an asymmetry factor up to A u ~ 10~ 4 B u , where B u = _B/10 14 G. In our case of the strange 
star in the CFL phase with a strong magnetic field the direct QDU process and formation and breaking 
of quark pairs are relevant as well as the reaction vv — > ee + near the star surface. Here we assume that 
the QDU processes for T < T ( ^f 6 i within the beaming angle (in accordance with the no go theorem 
O [46J) produce an asymmetry factor of the same order of magnitude as has been estimated in the 
literature. This A v factor then leads to a net momentum transfer from the neutrino flux to the star of 
mass M resulting in a time-dependent kick velocity 

V (t) = 4z r CFL dt ' L ^) > ( io ) 

which saturates at the magnetic field dependent asymptotic value 

v ~ 1.7 ■ B^M/MqY 1 km s~\ for A v ~ 10~ 4 B u , (11) 

as soon as the beaming ceases. Numerical results are shown in Fig. El for a typical strange star with 
M = 1.4 Mq and R = 10 km. In order to obtain large star kicks with v ~ 10 3 km s _1 , magnetic fields 
Bin, s ~ 10 17 G and large gaps A ~ 200 MeV are required. 
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Figure 2: Time evolution of the 
pulsar kick velocity for two val- 
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MeV (left panel) and A = 200 
MeV (right panel), for different 
magnetic fields B = B ins . 



4 Conclusions 

We have shown that in the presence of a strong magnetic field Bi nyS ~ 10 15 — 10 17 G, the early cooling 
evolution of strange stars with a color superconducting quark matter core in the CFL phase right after 
a short transport era is characterized by anisotropic neutrino emission, collimated within a beaming 
angle 9 U around the magnetic axis. Initial temperatures of the order of T ~ 30 — 60 MeV allow for the 
energy release ~ 10 53 erg. An energy ~ 10 52 erg is released within a narrow beam on time scales from 
10~ 2 s to 10 s, appropriate for the phenomenology of both short and long GRB. Exploring the beaming 
mechanism one may obtain a wide range of strange star kick velocities up to 10 3 km/s, in dependence 
on the magnetic field, pairing gap, radius, mass and burst duration of the star. The range of strange 
star kick velocities is thereby in a qualitative agreement with recent observational data on pulsar kick 
velocity distribution. All our estimates are very rough and essentially vary with the parameters of the 
model. Thus we just have shown a principal possibility of the application of the model to the GRB and 
kick description. A more systematic investigation, together with a consistent modeling of the compact 
star structure will be developed along the lines described in this contribution |37j . 

Acknowledgements 

J.B. and D.B. acknowledge support by DAAD grant No. DE/04/27956 and by the DFG Graduate 
School 567 "Strongly Correlated Many-Particle Systems" at University of Rostock. The work of D.N.V. 
has been supported by DFG grant 436 RUS 113/558/0-3 and by an RF President grant NS-5898.2003.02. 
H.G has been supported by DFG under grant 436 ARM 17/4/05. D.B. is grateful for the invitation to 
participate at the School for Nuclear Physics in Erice and for the support by DFG. 



References 

[1] I. Bombaci, Lect. Notes Phys. 578 (2001) 253 

[2] T. Piran, Phys. Rep. 314 (1999) 575 ; Phys. Rep. 333 (2000) 529 

[3] D.A. Frayl et al, Astrophys. J. 562 (2001) L55 

[4] G.E. Brown, C.H. Lee, R.A.M.J. Wijers, H.K. Lee, G. Israelian, and H.A. Bethe, New Astronomy 
5 (2000) 191 

[5] S. Rosswog, and E. Ramirez-Ruiz, AIP Conf. Proc. 727 (2004) 361 ; |astro-ph/030920Tl 



7 



[6] P. Haensel, B. Paczynski and P. Amsterdamski, Astrophys. J. 375 (1991) 209 

[7] V. Usov, ASP Conf. Series 190 (1999) 153 

[8] J.R. Salmonson, and J.R. Wilson, Astrophys. J. 578 (2002) 310 

[9] J.R. Salmonson, and J.R. Wilson, Astrophys. J. 517 (1999) 859 
[10] M.A. Ruderman, L. Tao, and W. Kluzniak, Astrophys. J. 542 (2000) 243 
[11] I. Bombaci, Nucl. Phys. A 754 (2005) 335 
[12] D. Aguilera et al, NATO Science Series 197 (2005) 377 

[13] G. Lugones, C.R. Ghezzi, E. de Gouveia Dal Pino, and J. Horvath, Astrophys. J. 581 (2002) L101 
[14] A.G. Lyne, and D.R. Lorimer, Nature 369 (1994) 127 
[15] H. Spruit and E. S. Phinney, Nature 393 (1998) 139 

[16] G. Hobbs, D. R. Lorimer, A. G. Lyne, and M. Kramer, MNRAS 360 (2005) 963 

[17] D. Lai, D.F. Chernoff, and J.M. Cordes, Astrophys. J. 549 (2001) 1111 

[18] M. Colpi and I. Wasserman, Astrophys. J. 581 (2002) 1271 

[19] V.S. Imshennik and O.G. Ryazhskaya, Astron. Lett. 29 (2004) 831 

[20] S.G. Moiseenko, G.S. Bisnovatyi-Kogan and N.V. Ardeljan, Proc. of IAU Coll. 192 (2003) 
[21] N.V. Ardeljan, G.S. Bisnovatyi-Kogan, K.V. Kosmachevskii, and S.G. Moiseenko, Astrofizika 47 
(2004) 37 

[22] Z. Arzoumanian, D.F. Chernoff, and J. M. Cordes, Astrophys. J. 568 (2002) 289 

[23] K. Rajagopal and F. Wilczek, At the Frontier of Particle Physics 3 (2001) 2061 

[24] K. Iida, T. Matsuura, M. Tachibana, and T. Hatsuda, Phys. Rev. Lett. 93 (2004) 132001 

[25] L. Lindbom and B.J. Owen, Phys. Rev. D 65 (2002) 063006 

[26] D. Blaschke, D. M. Sedrakian and K. M. Shahabasian, Astron. & Astrophys. 350 (1999) L47 

[27] D. Blaschke, T. Klahn and D. N. Voskresensky, Astrophys. J. 533 (2000) 406 

[28] D. Sedrakian, D. Blaschke, K.M. Shahabasyan and D.N. Voskresensky, Astrofizika 44 (2001) 443 

[29] I. Giannakis and H.-C. Ren. Nucl. Phys. B 669 (2003) 462 

[30] G. Baym, C. Pethick, D. Pines, and M. Ruderman, Nature 224 (1969) 872 

[31] E.M. Lifshiz and L.P. Pitaevsky, Statistical Physics, IX, Akademie Verlag Berlin, (1980) 

[32] D.N. Voskresensky, M. Yasuhira, and T. Tatsumi, Phys. Lett. B 541 (2002) 93 

[33] D.N. Voskresensky, M. Yasuhira, and T. Tatsumi, Nucl. Phys. A 723 (2003) 291 

[34] D. Blaschke, S. Fredriksson, H. Grigorian, A. Oztas and F. Sandin, Phys. Rev. D 72 (2005) 065020 

[35] LA. Shovkovy and P.J. Ellis, arXiv: hep-ph/0303073 (2003) 

[36] N. Iwamoto, Phys. Rev. Lett. 44 (1980) 1637 ; Ann. Phys. 141 (1982) 1 

[37] D. Blaschke, H. Grigorian, and D.N. Voskresensky, Astron. & Astrophys. 368 (2001) 561 

[38] D. Bandyopadhyay, S. Chakrabarty, P. Dey and S. Pal, Phys. Rev. D 58 (1998) 121301 

[39] Y.M. Loskutov and K.V. Parfenov, Vestn. Mosk. Univ., Fiz. 30 (1989) 3 

[40] A.B. Migdal, E.E. Saperstein, M.A. Troitsky, and D.N. Voskresensky, Phys. Rep. 192 (1990) 179 
[41] K.V. Parfenov, Sov. J. Nucl. Phys. 48 (1988) 651 
[42] K.V. Parfenov, Sov. J. Nucl. Phys. 49 (1989) 1126 
[43] A. Vilenkin, Astrophys. J. 451 (1995) 700 

[44] C. J. Horowitz and J. Piekarewicz, Nucl. Phys. A 640 (1998) 281 

[45] A.V. Kuznetsov and N.V. Mikheev, Phys. Lett. B 394 (1997) 123 

[46] A. Kusenko, G. Segre and A. Vilenkin, Phys. Lett. B 437 (1998) 359 

[47] J. Berdermann, D. Blaschke, H. Grigorian and D.N. Voskresensky, in preparation 



8 



